APS/ChauncerPRL211011 



Uncertainty Relation for Chaos 

A. YahalomQ M. Lewkowiczlll J. Levitanlil G. Elgressyi and L.P. Horwitfl 

Anel Umversity Center of Samaria, Ariel 40700, Israel 



o 

(N 
o 

Q 

On 



Ph. 
I 



> 
o 

(N 



X 



Y. Ben-Zion 

Department of Physics, Bar Ran University, Ramat Gan 52900, Israel 
(Dated: December 12, 2011) 

A necessary condition for the emergence of chaos is given. It is well known that the emergence 
of chaos requires a positive exponent which entails diverging trajectories. Here we show that this is 
not enough. An additional necessary condition for the emergence of chaos in the region where the 
trajectory of the system goes through, is that the product of the maximal positive exponent times 
the duration in which the system configuration point stays in the unstable region should exceed 
unity. We give a theoretical analysis justifying this result and a few examples. 
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There are many systems for which local instability 
leads to chaotic behavior. It is often possible to char- 
acterize this local instability in terms of the divergence 
of orbits detected by computing Lyapunov exponents [ij , 
or by geometrical methods in cases where a geometric 
description becomes available, such as in the application 
of the Jacobi metric and the time dependence of the re- 
sulting Jacobi equation for geodesic deviation [J-Q for 
Hamiltonian systems, or the criteria based on the cur- 
vature of the dynamical manifold obtained from the re- 
cently developed geometric embedding method (GEM) 
. These methods generally involve the computation of 
an exponential divergence as an indication of local in- 
stability. A given restricted region of exponential diver- 
gence, however, may not result in significant deviation of 
orbits, and therefore in chaotic behavior of the system. 
In this paper we give a quantitative bound which has the 
form of an uncertainty relation characterizing the effec- 
tiveness of the divergence in a locally unstable region. 
We find a relation between the time of passage of the 
orbit through an unstable region and the measure of in- 
stability, for example, a negative eigenvalue in the Jacobi 
equation, or the positive value of a Lyapunov exponent, 
which we denote in all these cases by A, which results in 
stability. This relation is of the form of a product 
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> 1, 



(1) 



similar to that of an uncertainty relation between time 
and frequency in optics, or in quantum theory. The basis 
for this relation is that if there is an indication of expo- 
nential divergence, it is the time the system spends in 
the region of instability which determines whether the 
exponential divergence has reached a significant magni- 
tude during this period. In general, the coefficient A is 
time dependent, but for a small region of instability, the 
relation equation ([T]) provides a useful measure. More- 
over, this is a necessary condition of instability, and if it 
is not satisfied by the trajectories under study, one can- 



not expect chaotic behavior. Formula ([1} is applied in 
several examples. 

In a recent study Q of the restricted three body prob- 
lem, involving several configurations of an Earth, Sun, 
Jupiter type system, we successfully applied the geomet- 
ric criterion of Q based on a geometrical embedding of 
the Hamiltonian orbits derived from a conformal trans- 
formation of the original Hamiltonian of the form 



H = 



P 

2me 



V{r,rj,ri2) 



(2) 



where 
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V{r,rj,ri2) = 

A ■> ■> r rj ri2 

(3) 

This potential has an adiabatic time dependence due to 
the motion of the large planet (Jupiter) through the pe- 
riodic dependence (frequency ujj of its radial variable Vj ; 
r is the radial coordinate of the Earth from the Sun, and 
ri2 is the Earth- Jupiter distance). Although the GEM 
criterion was developed for time-independent potential 
systems, it was shown in [s!] that the time dependence in 
this problem was negligible, provided the orbits did not 
pass too close to the boundaries of the physical region. 
Our prediction was that of stability with the exception 
of some small regions of instability near the foci (ap- 
sides) in highly eccentric orbits. These small regions of 
instability do not affect the overall stability of the sys- 
tem. In contrast, both Lyapunov methods and methods 
based on the time dependence of frequencies appearing 
in the application of the Jacobi metric, which predicted 
chaotic behavior for the system [3, indicated instability 
over large regions of the orbits. One therefore must un- 
derstand why exposure of the orbits to small regions of 
instability does not affect the outcome of the simulations, 
indicating completely stable behavior. These results were 
also obtained for the case of the Jupiter mass going to 
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zero, for which the GEM method predicts the stabihty of 
the two-body Kepler problem correctly. 

We have therefore investigated the application of the 
relation ([T]) , and found that it was indeed fulfilled in these 
cases. We have furthermore applied the idea to other 
known potentially chaotic systems in dynamical regions 
(according to the choice of parameters) where the passage 
time multiplied by the maximum negative eigenvalue is 
small (but not zero) and found that these systems re- 
mained stable in this regime. Moreover, we found in these 
examples that as the parameters are changed to induce 
chaotic behavior, the degree of chaos was well correlated 
with the growth of the uncertainty product. This is con- 
sistent with the interpretation that the passage time for 
exponential deviations in the orbits induced by the time 
dependent eigenvalues of the Jacobi equation are critical 
in the development of chaos for such Hamiltonian sys- 
tems. We have furthermore verified that a similar phe- 
nomenon occurs for the standard Lyapunov analysis in 
problems for which it is applicable. In the following we 
introduce a toy model for this phenomenon, and present 
results from simulations of motion based on a polyno- 
mial potential for the Hamiltonian system describing the 
Toda problem. We conclude with a discussion of the Ke- 
pler problem. 

Assuming that we have two trajectories with a differ- 
ence ^i, i — 1, 2, 3, ...n in n dimensions, we study a Jacobi 
equation with the simplifying assumption that we are 
close enough to the tangent space such that the covariant 
derivatives are well approximated by ordinary derivatives 
i.e., the geodesic deviation (Jacobi) equation is 



In the Lyapunov analysis: 
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(4) 
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in which Xi are the spatial degrees of freedom of the me- 
chanical system, V is the potential of the mechanical sys- 
tem and rrii are the masses of the particles. In the GEM 
analysis Q: 



r 1 

nH = — 



dV dV 



2{E — Y) dxi dxj dxidxj 



(6) 



For a slowly varying matrix N with a positive eigenvalue, 
^ may appear to be exponentially divergent on a small 
time interval. If we follow the orbit for some time, the 
maximum positive eigenvalue may go to zero, and the 
system will enter a stable regime. If the eigenvalue tends 
to zero sufficiently rapidly, the exponential divergence of 
the orbits, characterized by ^, would not be adequate to 
lead to chaos. To see how this phenomenon can develop 
quantitatively, we compute the solution to the time de- 
pendent equation Q in what follows. 



Let us define ^ = rji] equation (j4]) can then be written 
C = MC, (7) 

in which 
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and 
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where / is a n x n unit matrix. Obviously the value 
of C cannot diverge faster than what is dictated by the 
maximal positive eigenvalue of M over a certain duration. 
Let us now look at a specific example. For a 2 x 2 (2D) 
N we have a 4 x 4 M. The eigenvalues A of are then 
determined by: 

- XTtN + detN ^0 (10) 
The eigenvalues fi of M are determined by 

/ - fi^TrN -f det iV = 0. (11) 
so that fi^ = X. Then, local stability is determined by 

A > : ^ real, unstable 

A < 0, /i imaginary, stable (12) 

It might be convenient to take 

N = al + bai + ca^ {ai is a Pauli matrix) (13) 

as a real symmetric example. Then, 

a + c b , ^^^^ 



a — c 
and the eigenvalues become: 



A± = a ± y/p + c^. 
Now let us assume that: 

6 = pcos0, c = psin( 



(15) 



(16) 



where p > and 9 are some arbitrary parameters. In 
this form the eigenvalues will take the form: 



A± = a ± p. 



(17) 



Notice that the value of 9 does not effect the eigenvalues 
and hence we will take it as = 0. Moreover by proper 
temporal scaling we can always take p = 1. For a = 
— \a\ < 0, the eigenvalue A_ is negative and hence always 
stable; we thus need to analyze only the behavior of X^. 
The following form for a is assumed: 



t 

At 



(18) 
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FIG. 1. The case At = 5, the figure depicts £,i{t) (5i(0) = 1) 
by a red fine, £,2{t) (^2(0) = 0) by a blue hne while the black 
line describes an exponential growth according to the maximal 
exponent \+max ~ 1- (color online) 




FIG. 2. The case At = 0.1, the figure depicts ^i{t) (^i(O) = 1) 
by a red line, £,2{t) (^2(0) = 0) by a blue line while the black 
line describes an exponential growth according to the maximal 
exponent X+rnax = 1- (color online) 



for which At > is some given time interval. The eigen- 
values are then: 



Ah 
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At 



± 1. 



(19; 



In particular, A+ = 1 — ^ is positive for t < At (anc 
hence unstable) but negative for t > At and hence stable 
This model is therefore appropriate in order to study th( 
effect of the duration of the trajectory in the unstable 
region on the systems overall stability. The maximal in- 
stability exponent is \+max = 1 at t = 0. Thus equatior 
^ takes the simple form: 



At > 1 



(2o; 



as a necessary (but not sufficient condition) of instability 
To appreciate the significance of the above result let us 
look at two trajectories displaced at t = by an amount 
^ = (1,0) but with the same velocity rf = (0,0). In the 
case At = 5 depicted in figure [1] it is clearly seen that the 
exponential growth (black line) is not achieved by either 
component of ^. We see a growth by an order of magni- 
tude in the size of the displacement ^i, the displacement 
^2 achieves the same magnitude as although null at 
t = due to the coupling with ^i. The entire growth 
happens at an interval slightly longer than At = 5 in- 
cluding a slight overshot. Entering the stable region the 
displacements oscillate with a constant amplitude. Let 
us now consider the case in which equation (j20p is not 
satisfied, for example let At = 0.1; this case is depicted in 
figure [21 One can hardly notice the exponential growth 
in figure [21 as to and ^2 the fact that the "unsta- 
ble" region is so minute does not allow them to grow at 
all and they oscillate at an amplitude smaller than the 
original displacement. Having demonstrated the chaotic 
uncertainty principle with a toy model we now move on 
to more realistic examples. 

A generalization of the Toda potential is described by 
in this model the potential of the two-dimensional 
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FIG. 3. The product of A with the time spent in unstable 
regions At vs. energy. As the product becomes appreciable, 
the increase of chaotic behavior is seen clearly on Poincare 
plots, (color online) 



system is given by: 



1 



V{x,v) = |(a;2+j;2)^^2y 



7y' + ^2;4 + -y4 (21) 



Considering a unit mass, it was shown previously in Q 
that for energy E somewhat larger than 0.2 the system 
becomes chaotic. Calculating the value of AAt in which 
At is the duration spent in unstable regions, one arrives 
at figure [H It is clearly seen that the product of A with 
the time spent in unstable regions only reaches consider- 
able size when the system becomes chaotic, that is when 
E > 0.215. The Poincare plots for 0.2 < E < 0.215 does 
not show chaotic behavior since the product of A with 
the time spent in the unstable region is not consider- 
able. This is shown in figure [H The Lyapunov exponent 
just below the chaotic threshold oi E = 0.2 (not shown 
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Poincare plot for E=0.215 
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FIG. 4. The Poincare plot for E = 0.215. For this value 
XAt — 1.044 and a very mild chaotic behavior is seen as 
predicted by the chaotic uncertainty relation, (color online) 
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FIG. 5. A high eccentricity orbit, (color online) 



here) predicts chaotic behavior which is not manifested; 
this is another indication of the effectiveness of the GEM 
method. 

It is well known that the Kepler orbit which is clearly 
integrable and not chaotic, is predicted to be chaotic ac- 
cording to the Lyapunov stability criterion. At least one 
Lyapunov exponent is positive throughout the orbit in- 
dicating all over instability. In contrast, the geometric 
method developed by Horwitz et al. [6| generates the 
correct sign of the corresponding exponent. The same 
is true for the restricted three-body problem in which a 
"Jupiter" is constrained to move on a large circular tra- 
jectory Q. This system is described by equations ^ 
and ([3]). The above statement is true for low eccentric- 
ity orbits, however for high eccentricity orbits such as 
the one depicted in figure [5] even the GEM method suf- 
fers from local unstable exponents near the perihelion. 
Nevertheless, a calculation of the product of the expo- 
nent times the duration it takes the orbit to go through 
the unstable region indicates that the necessary condi- 
tion ([l} of instability is not met. For example the orbit 
of figure [5] exhibits a positive geometric exponent near 
the perihelion (being stable everywhere else) of the form 
described in figure [Bl For this unstable region we have 
Xmax^t c± 0.6, which is below what is required for chaotic 
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FIG. 6. The geometric eigenvalue of the Keplerian orbit near 
perihelion, (color online) 

instability. The Lyapunov criterion, however, implies in- 
stability over the entire orbit. 

Our analysis of the Kepler problem and the restricted 
three body problem using the GEM method has led us to 
investigate the meaning of localized instabilities. This in 
turn has led us to formulate the chaotic uncertainty prin- 
ciple of equation ([1} . The content of the idea was investi- 
gated through a toy model demonstrating that unless the 
condition is met one cannot expect small deviations to 
grow by more than an order of magnitude. The value of 
this idea was further underlined in the generalized Toda 
potential. Finally it proved most valuable in understand- 
ing the ineffectiveness of the very local indications of in- 
stability arising in the Kepler and three-body orbits that 
appear in the GEM analysis. We foresee the usefulness 
of the idea for studying chaotic systems in general. 
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